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Abstract 

Arnold pointed out that the Euler equation of incompressible ideal 
hydrodynamics describes geodesies on the group of volume-preserving dif- 
feomorphisms. A simple analogue is the Euler equation for a rigid body, 
which is the geodesic equation on the rotation group with respect to a 
metric determined by the moment of inertia. The metric on the group is 
left-invariant but not right-invariant. We will reduce the geometry of such 
groups (using techniques popularized by Milnor) to algebra on their tan- 
gent space. In particular, the curvature can be expressed as a biquadratic 
form on the Lie algebra. Arnold's result that motion of incompressible 
fluids has instabilities (due to the sectional curvature being negative) can 
be recovered more simply. Surprisingly, such an instability arises in rigid 
body mechanics as well: the metric on S0(3) corresponding to the mo- 
ment of inertia of a thin cylinder (coin) has negative sectional curvature 
in one tangent plane. 

Both ideal fluids and rigid bodies can be thought of as hamiltonian sys- 
tems with a quadratic hamiltonian, but whose Poisson brackets are those 
of a non-nilpotent Lie algebra. We will also describe a different point of 
view towards three dimensional incompressible flow in terms of the Cleb- 
sch parametrization. In this picture, the Poisson brackets are represented 
canonically. The hamiltonian is represented by a quartic function. 

This is meant mainly as an expository article, aimed at a mathemat- 
ical audience familiar with physics. Based on Lectures at the Chennai 
Mathematical Institute and the University of Connecticut. 



* rajeev@pas.rochester.edu 
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1 Introduction 



Unlike in mathematics, the problems of physics tend not to be very old. Exper- 
imental advances constantly invalidate old ideas or introduce completely new 
ones. However, the problem of understanding a non-integrable dynamical sys- 
tem ('chaos') is as old as physics itself and is still largely unsolved. A par- 
ticularly virulent example is the phenomenon of turbulence in fluid mechanics: 
when velocity exceeds a critical value, the flow suddenly becomes irregular and 
unpredictable except by very flne numerical methods. A theoretical understand- 
ing of this phenomenon, perhaps along the Hues of Wilson's theory of second 
order phase transitions, remains a great challenge. 

A separate question is whether the partial differential equations of hydrody- 
namics (Navier-Stokes or Euler) have unique solutions and how regular they are. 
This has been recognized as a mathematical challenge worthy of the best ana- 
lysts. Even though the two problems in physics and mathematics are different, 
one hopes that ideas from one will cross-fertiHze the other. 

We begin by reviewing the basic equations of the subject. For simpHcity and 
brevity, it is hard to beat the classic text by Landau and Lifshitz [1]. For the 
geometrical formulation, the references are the books by Arnold and Khesin [?] 
and Kambe [4] . The expository article by Milnor shows how the geometry of 
a group can be reduced to algebraic questions on its tangent space. This will be 
important for us because it allows us to avoid defining an infinite dimensional 
manifold. We can reduce everything to vector spaces and linear operators on 
them: much simpler technically, while the geometry of manifolds continues to 
provide powerful intuition in the infinite dimensional case. 



1.1 Ideal Fluids 

We will consider only ideal non-relativistic fiuids; that is, a fiuids whose veloci- 
ties everywhere are small compared to the velocity of light and in which losses 
due to friction (viscosity) are small enough to be ignored. 



1.1.1 The Two Time Derivatives 

There are two ways of thinking about the time dependence of any physical 
quantity in a fiuid: at a fixed location in space or, along the fiow of the fiuid. 
The first is the partial derivative ^ and the second is the total (or material) 
derivative ^. They are related by 

d_ _ d_ 

where v is the velocity of the fiuid and Cv is the Lie derivative. On a scalar 
field it is just 

d(t) d(j) , „ , 
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We will also be interested in the case of a density. A density on a manifold 
is simply a differential form of the highest possible rank ( i.e., equal to the 
dimension of the manifold) , pdx^ A ■ ■ ■ dx" . Although it has one independent 
component like a scalar, the Lie derivative of a density is different from that of 
a scalar: 

^vP = di[pv']. 
The divergence of a vector field can be defined as 

div V ~ -CvP ~ -dilpv'']. 
P P 

The time derivative of the velocity field itself requires a different idea. Since 
CyV ~ [w, v] ~ 0, the Lie derivative does not capture its variation due to the 
motion of the fiuid. Using a Riemannian metric g of the manifold in which the 
fiuid is moving ( more precisely, using its Levi-Civita connection V) we can find 
the acceleration of a fiuid element as ^ = + V„w 

dv^ dv 



dt dt 



where, as usual. 



1.1.2 Conservation of Mass 

The first law of motion of the fiuid is just the conservation of mass: 
We can write this also as 

dp 9p 

in terms of the material derivative. It is clear that this equation does not 
make use of the Riemannian metric of the manifold M in which the fiuid is 
moving. 

1.1.3 Conservation of Momentum 

For a fiuid without external forces, Newton's second law gives 

dV; 

where p is the pressure and ^ is the derivative taken in the co-moving 
reference frame of the fiuid. It includes the explicit time derivative as well as 
the change due to the motion of the fiuid element: 
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dv 1 

^ + v''VkV^ = V^p. (2) 

at p 
It can be re-expressed as the conservation of momentum density: 

d{pvi) 



dt 

where the stress tensor density is 



+ VfeTf = (3) 



=p5'i + pv'^v,. 

This equation does use the metric of the underlying space, which is often take 
to be EucHdean in physical appHcations. But there are physically interesting 
cases of fluid motion on a curved geometry as well: the ocean or the atmosphere 
of a planet can often be thought of as a fluid on the surface of a sphere. 



1.1.4 Equation of State 

So far we have one scalar equation and one vector equation for the unknown 
quantities p,p, ■y^.We need one more scalar equation to have enough information 
to determine them. This is given by an equation of state: a relation between 
pressure and density. A model that works well in many situations is the law 
p = Ap'^ , (polytrope) for some constants ^,7 characterestic of the fluid . For 
the atmosphere in the adiabatic approximation, 7 ~ 1.4. 



1.1.5 The Wave Equation of Sound 

Before studying any non-linear equation in depth, we must understand its linear 
approximation. If the gradient of the velocity is small and the departure pi of 
the density from some average value po is small, the above equations linearize 
to 



dp 
~dt 

diPQv') 
dt 



kV>i = 0. 



where 



K 



dp 



. By differentiating the first equation w.r.t. time and 

P=PO 

putting in the second we get 

which is the equation for some wave propagating with velocity c = y^- These 
are sound waves. Thus sound is the infinitesimal manifestation of fiuid fiow. 
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1.1.6 Isentropic Flows 

If there is a function w ('enthalpy') such that 

diW = -diP 
P 

the flow is said to be isentropic. It basically means that no heat is lost or gained 
by the system. Then the second equation of motion can be written as 

^ + v''^kV^ + ^'W = 0. (4) 

at 

We still need an equation of state giving w as a function of p. 
1.2 Incompressible Fluids 

If the time dependence of the density is small enough to be ignored, 

we say that it is incompressible. Then compressibHly k tends to infinity: the 
speed of sound is infinite. More precisely, the speed of the fluid flow is small 
compared to the speed of sound. In this case we get the equations of Euler: 

dk = (5) 



dt p 

The density and the metric of the manifold containing the fluid should be 
given; p is a Lagrange multiplier enforcing the condition of incompressibility [H 



Incompressible does not mean constant density 

It is often stated incorrectly that an incompressible fluid must have constant 
density. Incompressibility only means that the density is given as a function of 
space and is independent of time: it is not a dynamical variable. The atmo- 
sphere of the Earth, for example, has density decreasing with height; yet the 
atmospheric flows (winds) have velocities much lower than the speed of sound. 
The atmosphere is an incompressible fluid. 



1.2.1 Boundary Conditions 

We have not said much about the underlying manifold M on which the fluid 
flows. In most cases of interest it is a domain of Euclidean space. But it 
is useful to consider the more general case of a Riemannian manifold, as the 
mathematical concepts are more clear then. Even in physics, occasionally one is 
interested in the flow of a fluid on a curved manifold such as the sphere (ocean 
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currents). In general the manifold M will have a boundary. For an ideal fluid, 
it is suflicient that the velocity fleld be tangential to the boundary. 

In the real world, dissipation (viscosity) cannot be ignored near the bound- 
ary, even if it is small away from the boundary. Thus, the physically correct 
boundary condition is that the velocity field must vanish at the boundary, not 
just its tangential component. We will assume this stronger condition in what 
follows. 

Thus, all vector fields will vanish at the boundary. All dilfeomorphisms 
reduce to the identity at the boundary. 

1.2.2 Two Dimensional Incompressible Fluids 

A particularly interesting example is the case of incompressible fluids in two 
dimensions. Geophysical applications include ocean currents and atmospheric 
flows which are approximately two dimensional: the atmosphere and the oceans 
have a depth small compared to the diameter of the Earth. As an example we 
consider the case of flow in the plane and with p = 1 : a part of the sphere that 
is small enough that the curvature can be ignored. We will return to a more 
general theory later. 

We can eliminate w by taking the curl of the above equation. In two dimen- 
sions, the curl of velocity ('vorticity') is a scalar 

W = dlV2 — 92^1. 

Since 

di(vjdjV2) — d2 {vjdjVi) = diVjdjV2 — d2VjdjVi + VjdjU (6) 

= 9lWl9lW2 + 9lW2<92U2 — (92Wl9lWl — 92W292W1 + VjdJ^ 

= uj{djVj) + Vjdju (8) 
= -^'jdjUJ (9) 

Moreover, every incompressible vector fleld is of the form 

vi = d2X, V2 = -dix 
for some 'stream function' x- Vorticity is then 

w = -dfx - dlx = Ax. 

where A is a positive laplacian. Given appropriate boundary conditions, this 
is an invertible operator. Thus we can regard the vorticity as the dynamical 
variable and the velocity potential as derived from it by solving the above elliptic 
differential equation. 

Also, 

VjdjUj = d2xdiuj - dixd2Uj = {%, u} 

which is the Poisson bracket on a two dimensional phase space. It is anti- 
symmetric and satisfles the Jacobi identity. Thus, two dimensional incompress- 
ible flow reduces to the pair of equations 
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— + {x,w} = 0, u} = Ax. 

The Green's function of the laplacian (with appropriate boundary condi- 
tions), can be thought of a linear operator which solves the Poisson equation 
X = Kuj. Then the equation of motion of an incompressible fluid in two dimen- 
sions becomes the non-linear integro-diflerential equation 

^ + {KujM = ^- (10) 



2 The Rigid Body 

It is useful to start with an example of a mechanical system that looks like the 
opposite extreme from a fluid: a rigid body. We will see that there are many 
similarities in the basic mathematical formulation, although fluid mechanics 
is much more compHcated. It is interesting that the basic equations of both 
extremes are due to Euler. 



2.1 Euler Equations 

Recall the Euler equations of a rigid body on which no external forces are 
acting are an expression of the conservation of angular momentum L. In the 
non-inertial reference attached to the body itself, this takes the form 

— + f^xL = 

dt 

where Q, is the angular velocity. These quantities are related by a positive 
symmetric tensor (linear operator) , / the moment of inertia: 



L = m. 

The moment of inertia / is deflned in terms of the density of the rigid body 
as follows: 



It is obvious that the matrix M > 0. Suppose its eigenvalues are Afi, Af2, Afa, 
all positive numbers. Then M^k = Mi + M2 + M3 . The matrices M and / are 
diagonal in the same basis, and the eigenvalues of lij are 

h = Ah + A/3, h = A/i + A/3 I3 = Ml + A/2. 

Thus / > as well. The basis in which / is diagonal forms the principal 
axes of the body and /i, /2, /3are called the principal moments of inertia. 

If we compare with the Euler equations in two dimensions, we see that the 
vorticity is analogous to angular momentum and the velocity potential analogous 
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to angular velocity. Also, the elliptic differential operator A above is analogous 
to the moment of inertia. The cross product of vectors (which is anti-symmetric 
and satisfies the Jacobi identity) corresponds to the Poisson bracket of functions. 

L^o^, n^x, /^A,x^{,}. 

Of course, the angular momentum has only three independent components while 
vorticity belongs to an infinite dimensional space. So fiuid mechanics is much 
more complicated. 

There is a co-ordinate system, which moves together with the body, in which 
lij is diagonal: 

In terms of these the Euler equations become 



—1 + ^— ^{12^3^0, + ^—^Ogrii = —± + ^—^nin2^o 

at Ii at 12 at I3 

(11) 

They are solved in terms of the three Jacobi elliptic functions . 



2.1.1 Hamiltonian Formalism 

A rigid body moves so that the distance between any two points on it remains 
fixed. Thus its configuration space is the isometry group of R^. More precisely 
the connected component of the isometry group, which consists of rotations 
and translations. The translations are uninteresting as they simply describe a 
straight-line on which the center of mass moves. Thus, we can think of the 
group of rotations SO{3) as the configuration space of a rigid body. 

The components of angular momentum satisfy the Poisson bracket relations 

{Li,L2} = L2, {L2, -^3} = -^3, {L3,Li}~Li. 

Just as the kinetic energy due to translational motion of a particle is 
the rotational kinetic energy of a rigid body is given by 

t2 t2 t2 

(12) 

2/1 2/2 2/3 ^ ^ 

in the basis where lij is diagonal. Indeed, we can check that Euler equations 
above are implied by this hamiltonian and the above Poisson brackets: 
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2.1.2 The Elliptic Curve 



The Poisson algebra of angular momentum has a non-trivial center. That is, 
there is a polynomial in the generators that commutes with (has zero Poisson 
brackets with) all the generators: 

L^ = LI+LI + lI (13) 

It is the square of the magnitude of anguar momentum. In particular, it 
commutes with the hamiltonian and hence is a conserved quantity: 



Of course, the hamiltonian is itself a conserved quantity. 

Thus the solution to the Euler equations describes parametrically the curve 
which is the intersection of the sphere ifTSj) with the ellipsoid l|12p . This in- 
tersection can be either a union of disjoint circles immersed in or a single 
connected closed curve, depending on the values of H and L. 

So the solution of Euler's equation has a purely algebraic description. When 
solving algebraic equations, it is useful to continue to complex values, even if the 
physical values are real as in our case. The intersection curve is then a complex 
curve, a manifold of two real dimension. A moment's thought will convince you 
that this manifold must be a torus: it is compact, and the time evolution defines 
an everywhere non-zero vector field on it. (Only connectedness needs a proof 
which we skip.) Thus was born the theory of eUiptic curves. 

The theory of elliptic curves have been honed into a fine marble sculpture in 
the garden of mathematics. But some of the life has been lost in this process: 
the physics seems to be lost. 



2.1.3 Geodesies on SO{3) 

Euler equations have a natural geometric interpretation as geodesies on the 
rotation group, with respect to a metric determined by the moment of inertia. 
This metric is, in interesting cases, not bi-invariant. Instead, it is only invariant 
under the left action. Thus we can visualize the universal cover of the rotation 
group as a three dimensional elhpsoid. 

To see where this metric comes from, remember that the Lie algebra of a 
group can be thought of as the space of left-invariant vector fields. A positive 
quadratic form on the Lie algebra defines a left invariant metric on the group. 
Euler's equations describe how the tangent vectors to the geodesies evolve, with 
'time' having the meaning of arc length. 

This point of view will be very useful for us in understanding hydrodynamics. 
So we will describe in more detail the geometry of left-invariant metrics in a later 
section. 



9 



3 Hamiltonian Systems From Lie Algebras 



We will see that the Euler equation of both ideal fluid mechanics and the rigid 
body are special cases of a class of dynamical systems obtained from a real Lie 
algebra Q. We digress a bit to describe this class of systems. 

The set of observables of a classical mechanical system is a Poisson alge- 
bra. That is, a commutative algebra on which is deflned in addition a bilinear 
(Poisson bracket) such that 

1- {/i,/2} = -{/2,/i} anti-symmetry 

2. {{/i, M , h} + {{/2, h} , /i} + {{/3, /i} , h} = Jacobi identity 

3. {/i, hh] = {/i, /2} /3 + /2 {/i, /a} Leibnitz identity 

The algebra of functions f : Q* ^ R on the dual of a Lie algebra is an example 
of a Poisson algebra. The exterior derivative of such a function can be thought 
of as valued in the dual of Q* , which can be identified with Q. Thus it makes 
sense to take the Lie bracket of a pair of such exterior derivatives [d/i(a), d/2(a)] 
evaluated at some point a g C/*. A contraction with a itself gives a number. So 
we define the Poisson bracket to be 



{hj2]{a)^-ia{[dh{a),dh{a)]). 

The required properties follow from those of a Lie algebra and the exterior 
derivative. 

If we choose a particular function H as the hamiltonian, we get a dynamical 
system for which the time evolution of any observable is given by 

Now, every Lie algebra element lo £ Q defines a linear function on its dual. 
In this case, we get a simpler form of this equation: 

f + l<iH,.| = 0. 



3.1 Metric Lie Algebra 

A Metric Lie Algebra is a Lie algebra along with an inner product (usually not 
invariant). Now, an inner product is a symmetric tensor on Q, and hence can 
be thought of as a quadratic function on Q* . If we choose this function as the 
hamiltonian, dh : Q* ^ Q will be a linear function; which is another way of 
thinking of an inner product. 

If our Lie algebra Q admits an invariant inner product < , > (not necessarily 
H), we can write this equation a bit more explicitly. Then we can use the 
invariant inner product to identify Q* and Q and dH becomes just a linear 
operatoriC : Q* ^ Q : 
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H{uj) = i < uj,Kuj > . 
The equations of motion are then 

dcu 

An example of this is the rigid body, where the Lie algebra is the cross 
product, the invariant inner product is the dot product and K is the inverse of 
the moment of inertia, yielding Euler's equations . 

When the Lie algebra admits no invariant inner product it is not as easy 
to write an expHcit form for the equations. It is then useful to describe it in a 
basis Li of the Lie algebra. The Poisson bracket is determined completely by 
its effect on the generators: 

{L^,Lj} = c^jLfe 

where cf^ are the structure constants of the Lie algebra in this basis. Thus 
for example, 

{Li, LjLm} = c'^jLkLm + Ci„^LjLk 

We are thinking of the basis elements of G as generators of the algebra of 
polynomials on G*- This is a formal translation of the point of view of most 
physicists. 

A particularly simple polynomial is given by an inner product on the Lie 
algebra. Expressed in the above basis, 

H = -W^LiL, 

where K^^ = /i-'*. If we use this as a hamiltonian, the equations of time evolution 
can be written as 

If the inner product h where invariant, it would satisfy the identity 
and there would be no time evolution. 
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3.2 The Metric Lie Algebra of a Two Dimensional Incom- 
pressible Fluid 

The set of functions on the plane form a Lie algebra under the Poisson bracket 

{f,9} = d2fdig~ difd2g. 
There an invariant inner product in this Lie algebra 



< /,5 >= 




If choose as hamiltonian the positive quadratic form 

H{iLj) = i < cj,Gcj > 

where G is the Green's function of the Laplace operator A, the equations of 
motion obtained are exactly those of an incompressible fluid in two dimensions. 


Thus the inner product that defines an incompressible fluid is 
J f{x)G{x,y)g{y)dxdy 

where 

{df + di)G{x,y) = -S{x,y) 
with appropriate boundary conditions. 

3.3 Dimensions Greater Than Two 

Even if the dimension is greater than two, the hamiltonian of incompressible 
fluid is still a quadratic function on the Metric Lie Algebra of incompressible 
vector flelds. However there is no longer any obvious invariant inner product, 
so we have to contend with describing the equations a bit indirectly. 
The condition of incompressibility 

dk {pv") - (14) 

is equivalent to 

Cvp = 0. 

It follows that the commutator of two incompressible vector fields is also 
incompressible: the set of solutions Vp of[l4]is a Lie algebra with the same Lie 
bracket as before. A linear function on Vp is of the form / pv'^a.idx for some 
1-form a. But, a gauge transformation 

a 1-^ a + dA 

^We must assume appropriate boundary conditions on u> so that H{uj) exists. 
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leaves this function unchanged. Thus the dual of Vp is the space of 1-forms 
modulo exact 1-forms: 

v; = AVdA°. 

The kinetic energy of the fluid is a simple quadratic form on the Lie algebra 
of incompressible vecor flelds.: ^ / pv^v^ gijdx. The dual variable dual to can 
be thought of as just 

the square brackets being there to remind us that we must take the equiv- 
alence classes under gauge transformations. Thus the Hamiltonian can only 
depend on the vorticity 

uj = da, LUij — dittj — djai. 
But it must be a zeroth order operator in terms of velocity. 

where G is the Green's function of the elliptic system 

UJ = da, di (^pg^-'aj^ = 0. 
The Euler equations in vorticity form are 

where we are to regard v as the unique incompressible vector fleld determined 
by the vortcity by 

di [pv'] = 0, 9, [gjkv''] - dj [gikv''] = Uij 
with appropriate boundary conditions. 

4 Three Dimensional Incompressible Flow 

This is the most important kind fluid flow: the vast majority of physical phe- 
nomena take place in three dimensions and the velocities are small compared to 
the speed of sound. 

4.1 The Clebsch Variables 

In this case Euler equations can be expressed in a simpler form using a parametriza- 
tion due to Clebsch |^: 

(jj = dq A dp (15) 
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It is clear that this is only possible because dw = 0. Locally, any one-form in 
can be expressed in the form 

a = d\ + qdp 

so that any exact two form is locally of the form [l5l 

In terms of these variables, Euler equations become the statement that p, q 
are constant along streamlines: 

% + v'd.p = 0, § + v'd.q = 
ot at 

The velocity is determined in terms of p, q by 

^g^^[djX + qd,p] 
Here, A is ehminated by the condition of incompressibility 

d,[pv'] = ^ d,W^{dj\ + qd,p)] = 0. 
4.1.1 Canonical Relations 

The Clebsch variables p, q are canonical conjugates of each other. That is, if we 
postulate canonical commutation relations 

{p{x),q{y)} = 5{x - y), {p{x),p{y)] = = {q{x),q{y)} 

and the hamiltonian 



with wMetermined in terms of p, q as above, we get the Clebsch form of the 
Euler equations. 

To see this note that the canonical commutation relations give 

{F,p(2/)} = --^, {F,q{y)]= 



5q{yy ' • 5p{y) 
for any function of p, q. For example, if 

F = I rd.qdjpdx 



{F,p{y)} = -dj^d^p, {F, q{y)} = -dj^d^q (16) 

if /'^ is independent of p, q. Now, the kinetic energy of the fluid can also be 
writen as 

H J pv^v^g.jdx ^2 J ^''^^^3'^^ \ J f'^^^I^JP 
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where pv^ = djf^^ and = 



— is the velocity potential. Therefore 




which are the Euler equations. 

An extra factor of 2 appears (cancelling the ^ in the Hamiltonian) compared 
to [m because Z*-' itself depends linearly on v and hence on p, q. 

4.1.2 The Moment Map 

Marsden and Weinsetin [6] showed that the Clebsch parametrization has a nat- 
ural geometric interpretaion. On the space J- of real valued functions on a 
manifold with density, there is an inner product 



This can be used to turn J- O J- into a symplectic vector space: the conjugate 
pairs of functions p, q parametrize this phase space. Since volume preserving 
diffeomorphisms preserve the inner product above, they must act as canonical 
transformation on this phase space. The infinitesimal generaror these transfor- 
mations is just vorticity. The Clebsch parametrization 



in Euclidean space. The dot product is a rotation invariant inner product 
in R^, which turns © into a phase space on which the rotations act as 
canonical transformations generated by L. 

A direction for further research 

This point of view allows us to generalize the theory of a three dimensional in- 
compressible fluid to the case where the underlying manifold is non-commutative 
sphere: a kind of 'regularization' of hydrodynamics where each point is replaced 
by a fuzzy object which is the average of many points. In addition to being more 
mathematically rigorous (the equations of motion are finite dimensional ODEs 
instead of PDEs) this might be a physically realistic description of the large 
scale behavuor of fluids: the small scale fluctuations in velocity are averaged 
out to get an effective theory that s not no longer local. 




fgpdx. 



Lo = dq A dp 



is analogous to the formula for angular momentum 



L = 



r X p 
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4.2 The Lie Algebra of an Ideal Isentropic Fluid 

The set of vector fields V on a manifold is a Lie algebra under the commutator 
or Lie bracket, 

V acts on the space of scalars J- through the derivative 

The sum of the two vectors spaces V ® T ^ Q \s thus a Lie algebra as well: 

{u,(j)) 

This is the semi-direct sum of the Lie algebra of vector fields and the abehan 
Lie algebra of scalars. 

The dual of J- is the space of scalar densities J-* : 




The dual of V is the space of co- vector densities V* : 




Given functions F, G : V* ® JF* ^ i?, their Poisson brackets are given by 

[ fSF^SG SG ^ 5F\ ^ f.fSF^SG SG ^ 5F\ ^ 
J yojk op djk dp J J \djk oji djk oj^J 

This can be written also as 



SF 



SF 



5p 



SF 



Sji 



dx - I -^dk 



SF_ 



dx. 



Choosing G to be a linear function, we get a particularly convenient form of 
the Poisson bracket: 



5p Sji 



{F, p} = -dk 



SF_ 
^Sjk 



So far we have not needed any additional geometric structures such as a 
Riemannian metric: all the derivatives above make sense without the need of a 
connection. 
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Now we identify the variables ji , p physically as momentum density and 
mass density respectively. The Hamiltonian is the sum of kinetic and potential 
(internal) energies of the fluid: 



U{p)dx 



where is a Riemannian metric on the underlying space. 
On physical grounds we know also that momentum density is mass density 
times velocity: 



so that 



Ji = gikpv 



The equations of motion 

at 

become 



V . 



dp 



IgHj^ dU_ 

2 p2 + Qp 



We can simplify the first of these 

dji 1 

+ Vi {jkv') -pdk [gijv'-v^] - pv'gijS/kV^ - pdk 



dU_ 
dp 



where Vi is the covariant derivative of the Riemann metric g.y.The first two 
terms on the right hand side cancel each other. If we identify the enthalpy as 



w 



dp 



we get 



or 



'dt 



w'^VfcW* + g'''dkw = 



along with the conservation of mass: 
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These are exactly the Euler equations for an isentropic ideal fluid we obtained 
earlier, except that there we looked at the special case of the Euclidean metric. 
The Lie algebra V O J- as well as the Poisson algebra following from it are 
independent of the choice of metric. But the hamiltonian H depends on this 
choice. In this case, the hamiltonian is not a quadratic function, unHke in the 
case of the rigid body or the incompressible fluid. If the equation of state is a 
polytrope w = it is possible to think of the hamiltonian as a cubic function. 
Perhaps this is worth deeper study. 

5 Riemannian Geometry 

We review here some facts about Riemannian manifolds, to rephrase some basic 
facts in a language convenient for our purposes. This is not meant as an intro- 
duction to Riemannian geometry. There are several excellent texts available, in 
particular the one by Chavel [9]. 

5.1 Covariant Derivative 

A covariant derivative (connection) V„w of a vector fleld v on a manifold M 
along another vector fleld u should satisfy the conditions of linearity in u, v as 
well as 

V„[/u] = /V„u + VfuV = fVuV. 

Thus it involves flrst derivatives of v and no derivative of u. Explictly in 
co-ordinates 

for a set of connection coeflicients F^j,. We will only be interested in con- 
nections without torsion: 

WuV — V^w = [u, v\. 
The cuvature is the tensor fleld deflned by 

R{U,V)W = Vyu,v]W - VyVuW + VyVvW. 

Given a Riemannian metric g on the manifold, there is a unique connection 
of zero torsion and which preserves the metric: 

g{VuV, w) + g{v, Vuw) u {g{v, w)) . 
ExpHcitly in co-ordinates, this connection has coefficients 
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rj-fc = i^g" {djgik + dkgij - dig^k} 



The curves that extermize the action 



or 



g{i, x)dt 

are the geodesies; they satisfy the equation 

Vii = 



0. 



It will be convenient to define the curvature bi-quadratic form 

R{u, v) — g {R{u, v)v, u) 

5.2 Geodesic Deviation and Curvature 

The curvature form determines the rate of deviation of nearby geodesies from 
each other. 

Consider a geodesic x^{t) , expressed in a co-ordinate system centered at the 
initial point: x*(0) = 0. 



j dx^ dx 

^^^~dr 



0. 



Dx 
~dt 



0. 



Let the initial velocity be w G TqM. An infinitesimally close geodesic to this 
one will satisfy 



2„,i 



dx^ dx 



dt^ ' ^'^''■^ dt dt ^ 
D^y 



0. 



df^ 



+ R{x,y)x = 0. 



The vector field y along the original geodesic connects the points at equal 
time on two nearby geodesies. It is called the Jacobi vector field. We can derive 
an equation for the length squared of the Jacobi field 



2 dt^ 



Dy 



dt 



R[x,y). 



Suppose the initial conditions for the Jacobi equation are 
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y^^^ = ^^ = * = 

That is, we consider two geodesies starting at the same point but with 
slightly different initial velocities. 
Then 

\y{t)\^=t'\u\^-^^R{u,v) + 0{t% 

The first term would have been the answer in Euclidean space. Thus, if 
R{u,v) < 0, the geodesic with tangent vector v is unstable with respect to an 
infinitesimal perturbation in the direction u. 



5.3 Curvature as a Biquadratic 

We saw that the curvature tensor on a Riemannian manifold describes the be- 
havior of geodesies under small changes of the initial conditions. Therefore, 
it controls the stability properties of a physical system whose time evolution 
is along geodesies. It will be useful to calculate this tensor for left-invariant 
metrics on a Lie group to understand the stability of systems such as the rigid 
body or an ideal fluid. The formulas can get quite complicated, but a trick 
mentioned in Milnor's article allows us a simpler description. We will derive a 
simple formula that is quite useful in our appHcations by going a little beyond 
Milnor in this direction. 

Recall that a co-variant symmetric tensor on a vector space is exactly the 
same thing as a quadratic form. A covariant symmetric tensor is a bilinear map 
Q : V X V R, satisfying Q{u,v) = S{v,u). A quadratic form is a function 
Q : V ^ R that satisfles the scaHng property 

Q{Xu) = \^Q{u). 

Given a covariant symmetric tensor we can construct a quadratic by taking 
the special case when its entries are equal: 



Q{u) = Q(m, u). 

Conversely, we can get a symmetric tensor from a quadratic by 'polarization': 

^, . Q{u + v)-Q{u)-Q{v) 
Q[u,v) = . 

These maps are inverses of each other. 

In the same spirit, a tensor with the symmetries of the Riemmanian curva- 
ture is fully determined by a bi-quadratic form on the tangent space. In the 
appHcations we have in mind, this is a much more convenient description, as 
we will be able to write expicit formulas and identify their positivity properties 
more easily. 
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Recall that a bi-quadratic is a function on a vector space T : V x V ^ R 
that satisfies 

T(Au, v) = )?T{u, v), T{u, Xv) = X'^T{u, v). 
Now, the Riemann curvature tensor is a fourth rank tensor defined by 

r{u, V, w,x) = g {V[u,v]W - V„V„w + V„V„w, x) . 
It defines a biquadratic on the tangent space by choosing w = u,v = x [8] 

r{u, v) = g (V[„,„]U, v) - G (V„V^„u, v) + G (V„V„u, v) 
This function satisfies 

r{Xu,v) = X'^r{u,v), r{u,v) = r{v,u), r{u.u) = 0. (17) 

The meaning of these conditions is that k(u, v) = —; — ^ r(u,v) — depends 
on the subspace defined by u, v. To prove this, we just have to show that the 
condition 

r{au + bv, cu + dv) = [ad — hcf'r{u, v) (18) 

is equivalent to the conditions ifT7|l . The denominator is the area of the par- 
allelogram in the tangent space with sides u, v. The ratio fc(u, v) is the sectional 
curvature of the plane spanned by w, v. 

As Milnor points out, the conditions (fT7| imply in turn all the symmetry 
properties of a Riemann tensor. Also, k{u, v) determines fully the Riemann 
curvature tensor through a 'polarization' identity [9]: 



r(w, V, w, x) 



6 



dsdt 



{r{u + sw, V + tx) — r{u + sx, v + tw)} 



s=t=0 



The Ricci tensor is the trace of the Riemann tensor; equivalent to it is the 
quadratic form 

r{u) = r(u, Ci, u, ej)g^-' 

where e^is some basis in the tangent space and 5'-' is the inverse of the matrix 
of inner products gij = g{ei, Cj). 

From our current point of view it can be viewed as the average over all 
vectors 



r{u) — J r(u,v)d^ig{v). 

The average is with respect to the Gaussian measure on the tangent space, 
with zero mean and covariance given by the metric g. Recall that in terms of 
components 
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Also, the Ricci scalar is the further average over u: 




This point of view has the advantage that it extends to infinite dimensions. 
Also it is more natural in applications involving random forces. 

6 Geometry of Left-Invariant Metrics 

Let be a Lie group. Its Lie algebra G can be thought of either as the tangent 
space at the identity or as the space of left invariant vector fields on ®. An inner 
product G on is thus equivalent to a left-invariant Riemannian metric on 6. 
We will study the geodesies and curvature of this metric. Using homogenity, 
all computations can be reduced to the Lie algebra. The basic reference is the 
article by Milnor [8], especially Section 5. We will go beyond Milnor in deriving 
explicit formulas in a form useful for applications to mechanics. 

6.1 Covariant Derivative 

The covariant derivative (Levi-Civita connection) is determined by the deriva- 
tive of a left-invariant vector field by another. We will denote the covariant 
derivative on a group manifold by D to distinguish it from the covariant deriva- 
tive on a general Riemannian manifold, which we denote by V. This will be 
helpful when we talk of the diffeomorphism group of a Riemannin manifold: the 
derivative on the group of diffeomorphisms and that on the underlying manifold 
are closely related, but not identical notions. 

The conditions of zero torsion and preserving the metric become 

DuV - DyU = [u,v\, 
G{D^v,w)+G{v,D^w)=0; 
The zero torsion condition can also be written as 

G {DuV, w) — G {DyU, w) — G {[u, v],w) . 
Its cyclic permutations give 

G {DyW, u) — G {D^v, u) — G {[v, w], u) 

G [D^u, v) ~ G {DuW, v) ^ G{[w, u\,v) 

Adding the first and third and subtracting the second, and using the invari- 
ance of the metric 
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GiD^v,w) ^^{G{[u,v],w)~Gi[v,wlu) + Gi[w,u],v)} 
Define the linear operator u : G ^ G hy 

G{uv,w) ^ G{[u,vlw) + G{v,[u,w]) . (19) 

If the metric G were invariant under the Lie algebra, u would vanish for all 
u. Thus it measures the deformation of the metric. Now, 

G {uv, w)+G {vu, w) = G {[u, w],v) + G{[v, w],u) . (20) 

Thus 

DuV = ^{[u,v] - uv ~ vu} . 

6.2 Geodesies on a Group Manifold 

Given a curve 7 : [a, 6] ^ (5, its tangent vector at each point can be thought of 
as an element of the Lie algebra: 

d"/ 

-dt^^""- 

Thus the tangent vectors give a curve in the Lie algebra v : [a, b] G- 
We can defines the action of the curve as 

S{l) = \f^G{v,v)dt. 

(Some geometers call this quantity the energy; physicists would call it the 
action.) 

A geodesic is a curve at which the action is stationary with respect to small 
variations. 

To compute this variation, let us consider a one parameter family of curves; 
that is a map 

4> : [0, e] X [a, 6] ©, for some positive number e. We require that the initial 
and final points are left unchanged 0(s, a) and (/)(s, b) are independent of s. 
Define 

So u(s, a) = u{s, b) = 0. 

We get the integrability condition 

dgV = dfU + [v, u]. 

Regarding (f>{., s)for each value of s as a curve, the action becomes a function 
of s. Its derivative is 
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dsS{js)= / G{dsv,v)dt= / G{dtu+[v,u],v)dt 

J a J a 

Using the linear operator v : Q —f Q defined earlier 
G {[v, u],v) = G (u, vv) . 

By integration by parts 

dsS{-fs) = G{u,dtv - vv)dt. 

J a 

Thus the geodesic equation on the group becomes the ODE on the Lie algebra 

dtv + D^v = 0. (21) 

6.3 Curvature of a Left-Invariant Metric 

We will now calculate explictly the curvature bi-quadratic for a Metric Lie Al- 
gebra; i.e., for a Lie algebra {G, G) with an inner product on it. 

G (DyDuU, v) — ~G{DuU, Dyv) — — G {uu, vv) . 
-G{D.aD,u,v) ^GiD,u,Duv) ^ ~\\[u,v]f + huv + vuf. 



G {D[^^y]U, v) = ^G {[[u,v],u],v) - ^\[u,v]\^ + ^G {[v,[u,v]], u) 



Thus 



3 1 1 

R{u,v) ~ —^\[u,v]f+- {G {[[u,v],u],v) + G{[v, [u,v]],u)}+-\uv+vuf—G{uu,vv) 



The middle term can be further simplified using 

G([u, w],v) = G{uw, v) — G(w, [u, v]) 
and G{uw,v) = G{uv,w). We get 



i {G{[[u,v],u],v) + G{[v, [u,v]],u)} = + ^G{[u,v],vu - uv) . 

Thus 



R{u, v) = -| [m, u]p + —G{[u,v],vu — uv) + — \uv + vu\'^ — G {uu, vv) 
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If the metric is bi-invariant, only the first term survives: the curvature of 
a bi-invariant metric is postive. We can 'complete the square' on the first two 
terms to get 

R{u, v) = - \[u,v] + vu — uv\'^ + G{uv, vu) — G {iiu, vv) (22) 
This simple formula appears to be a new result. 

We already see something important: the sectional curvature in any plane 
that contains a Killing vector is positive: 

u = => R{u,v) = -\[u,v] + vu\'^ . 

This result is known but is hard to see with the formulas known in the 
literature. 

6.4 Example: The Two Dimensional Lie Algebra 

The only non-abelian Lie algebra of dimension two has commutation relations 
in terms of basis vectors 

[eo,ei] = ei 

or in terms of components 

[u,v] = {O,uovi - vqUi) . 

The only metric is (up to a change of basis that does not change the above 
commutation relations) 

G{u, v) ~ uqVo + uivi. 
The corresponding group manifold is the half plane, 

H = {{xo,xi)\xo > 0}} 
with the multiplication law 

{xo,xi){x'q,x[) = {xoXq,xox[ +xi). 

The corresponding left-invariant metric is the Poincare metric. It is well 
known that this metric has constant negative curvature. Hence it can serve as 
a simple model of an unstable dynamical system. 

We get 

G{vv, w) = G{v, [v, w]) = viivQWi — wqvi) 

so that 

vv = (-Wi,UoUl). 
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The geodesic equation in the Lie algebra becomes 

dvo 2 d.vi 



dt dt 



The solution is (scaHng t so that \v\ = 1) 

1 



vq = — tanhi, vi = 



coshi 

To get the curve in the group we must solve 

dg 

where gv denotes the action of the group on its Lie algebra. 
Viewing v as an element infinitesimally close to the identity 

Thus 

dxo dxi 

_=^o.„,_...o.,. (24) 
The Riemannian metric is thus the Poincare' metric 

^^2 _ (ixQ + dx\ 

Xq 

Thus geodesies are 

( . ,f , y B + Atanh[t - to]). 
cosh(t — to) 

(It is convenient to choose the origin of time as the point where vq (t) vanishes 
instead of the starting time.) These are the well-known semicircles 

xlit) + {x,it)-Bf^A'. 

6.4.1 Conserved Quantities of Geodesic Motion 

The Killing vectors of the Poincare' metric 

fa ^ xodo + xidi, /i = 9i, f-i ^ (xf - xl)di + 2xoxido 
form the SC{2, R) Lie algebra: 

[/07/1] — fii [/o,/-i] = /-I, = 2/0- 

The corresponding conserved quantities of geodesic motion are 

Fa = xqPq + xipi, Fi=pi, F-i = {xl - xl)pi +2xoXipo 
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where Pa are the canonical conjugates of x° : 

{Pa.Xb} = Sab- 

Since the geodesic equations have as hamiltonian 
we get = g°'^Pb- Expressing in terms of the variables Waintroduced earlier 

Pa = — • 

Thus we see that the conserved quantities of the equations I|23l24p above are 

Xq Xq Xq 

Of these, Fq, Fi generate the left translations of the half-plane thought of as 
a group. Now, note that Fq — FiF-i = + v\ which is an 'obvious' conserved 
quantity. Putting in the formula for the L's we get the semi-circle for the shape 
of the geodesic. 

Curvature Form 

It is straightforward to calculate the curvature biquadratic for this case: 

R{u,v) = - [uqVi - VqUi]'^ 

which simply says that the sectional curvature is: 

K{u,v) = -1. 

We use this simple case to check the numerical factors in our formula for 
curvature. 

6.5 Geodesies on S0{3) 

Consider as a Lie algebra with the cross-product as the Lie bracket: 

[u, v] = {U2V3 — U3V2,U3Vl — U1V3, V2 ~ U2V1) . 

A corresponding Lie group is S'0(3). Any inner product in can be brought 
to the diagonal form by a rotation without changing the Lie bracket: 

G{u,v) = Giuivi + G2U2V2 + G3U3V3 

Thus 

G{uv, w) = Gi {[M2W3 — U3V2] wi + [U2W3 — U3W2] vi} 
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+ G2 {[usVi - U1V3] W2 + [U3W1 - U1W3] V2} 



+ G3 {[uiV2 - U2V1] Ws + [uiW2 - U2W1] Vs} 

SO that 

r~ 1 Gi — G3 , G2 — Gi 

Ml = — 7; "ai's H U3V2: ■ ■ ■ 

The dots denote three more relations obtained by cyclic permutations. In 
particular, 

r~ 1 G2 — G3 

Ml = — 7^ — V2V3 



The geodesic equations become 

dvi G3 - G2 



V2V3 = 0, ■ 



dt Gi 

Calculating as above gives the formula for curvature 



n/ N (G2 — G3) + 2Gi(G2 + G3) — 3Gi 2 

V) = — [U2V3 - V2U3) H 

4Crl 

In particular, if we choose u = (0, --^,0)to be a unit vector in the second 
principal direction and v to be a unit vector in the third direction the sectional 
curvature of the 23— plane is 

^ ^ (G2 - Ga)^ + 2Gi(G2 + G3) - 3Gf 
4G1G2G3 

The others are given by cyclic permutations. 
6.5.1 Stability of the Rigid Body 

Even for an anisotropic rigid body with Gi < G2 < G3, the rotations around 
the principal axes are time- independent: (vi,0,0) for example is a solution of 
the Euler equations [l]. Under small perturbations, the rotations around the 
principal axes with the largest and smallest moment of inertia are stable. But 
a rotation around the principal axis with the middle value of moment of inertia 
is unstable: small perturbations grow exponentially. 

The geometric pitcure allows us to generalize this analysis to time-dependent 
solutions. The geodesic deviation equation shows that a small perturbation 
along u to a geodesic with tangent vector v will grow exponentially if R{u, v) < 0. 
Let us see under what conditions this is possible. A change of variables from 
the principal moments of inertia to the principal curvatures of the body will 
help us. 
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Recall that the principal moments of inertia are 

Gi = 7\/2+A/3,••• 
where Afi, Af2, A/3 are the (always positive) eigenvalues of the moment ma- 
trix M = j px ® xdx. 

Define their reciprocals 

1 r /• ]-i 

= — = / p{x)xldx 

In terms of these, 

(G2 - G^f + 2Gi(G2 + G3) - icl = A[MiM2 + M^M^ - Msi-fal ■ 
4 

= [t^2 + A^s - Ml] 



R(U,V)^ -, ■ ^[U2V3~U3V2\ H 

mi(m2 + m) 

Thus, if the three principals pi, p2, ps satisfy the triangle inequalities, 

Pl+P2>P3, M2 + M3>Ml, M3 + Ml>M2 

the curvature R{u, v) will be positive for all pairs: a rigid body that is not 
'too anisotropic' is geodesically stable. But these inequalities can be violated 
for a 'fiat enough' shape. Let us look at an example. 

Stability of a rotating cylinder 

If the body has an axis of symmetry (say the third axis) pi = p2 and 

R{u, v) = { [U2V3 - U3V2]'^ + [U1V3 - M3Ul]^} 

pi{pi + P3) 

H 7; Pl'^2 - U2V1] 

2piP3 

So an instability arises if 

P3 > 2pi. 

Consider a rotating cylinder of uniform density. Its height is h and its radius 
is r. It is elementary that 

Ml = A/2 — , M3 = — . 

4 12 



Ml = A'2 = — , P3 



r 



Thus a cylinder with too small a height 
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will be unstable. If it is rotating around one of its diameters, it is unstable 
with respect to a change of axis towards one of the other diameters. An example 
of such motion is a coin toss. 



7 The Geometry of DifFeomorphisms 

It is a remarkable fact that the set of Diffeomorphisms of a Riemannian manifold 
is itself a Riemannian manifold. This "higher Riemannian geometry" is the 
proper language of many interesting physical systems such as ideal fluids. Such 
repetitions of structures at a higher level happen quite often in mathematics: 
the set of Riemannian metrics on a manifold is itself a Riemmannian manifold; 
the set of complex structures is a complex manifold; the set of Kahler structures 
is itself a Kahler manifold and so on. 



The simplest manifold is the circle. So the first example of a Diffeomorphism 
group must be Dif f{S^). The standard metric on the circle leads to the L^— metric 
on this group. 



The deformation tensor of a vector field on the circle is easily found from 



G{uv, w) — G{uv'~vu' , w)+G{v, uw'—wu) — / {uv'w — vu'w + vuw' — vwu'} dx. 



7.1 The Circle 




to be 



u = ~3u' . 



Thus the covariant derivative on the Lie algebra is 



DuV = — {[u, v] — uv — vu} = u'v + 2m?;'. 
The geodesic equation is the inviscid Burger's equation 




The curvature form 



R{u,v) = - \[u,v] + vu — {twp + G(uv,vu) — G {uu,vv) 
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reduces to 



R{u,v) = \[u,v]\^. 



Thus the L^-metric on Diff{S^) has positive curvature. 
A Technical Remark 

This contradicts the statement in [4] and in several other places in the math- 
ematics literature that the L^— metric on the group of diffeomorphisms of M 
has zero curvature if {M,g) is itself flat. In particular, it is not possible to 
calculate the curvature of the subgroup of incompressible diffeomorphisms by 
using the Gauss-Codazzi formula as is often claimed. This is despite the fact 
the final answer has an uncanny resemblance to this formula . The confusion 
arises because there is a connection on Dif f[M) whose curvature form is the 
average of the curvature form of M. But it does not preserve the metric on the 
group manifold: 



where dvg is the Riemannian volume element on M . We take a more direct 
route in the next section to avoid this pitfall. 

7.2 Incompressible Diffeomorphisms 

Given a scalar density on a manifold, the set of diflteormorphisms that preserve 
it form a subgroup: 



Its Lie algebra is the set of vector fields of zero divergence. Recall that the 
divergence of a vector field is defined with respect to a density as 



it follows that incompressible (divergenceless) vector fields form a sub-Lie 
algebra, which we will call Q. 

Given a Riemannian metric g on M there is an inner product on the space 
of vector fields on M: 




Dp = {0 e S 1/9(2:) = detdcj){x)p{(f){x))} . 



divu ^ —di [pu'l . 
P 



From the identity 



div[u, v] 



u [divw] — V [divu] 
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We will call this the L^— inner product even though we are not completing 
the space of vector fields with respoect to the induced norm. Define now the 
deformation tensor u of a vector field: 

G{uv,w) ^ / {g{[u,v],w) + g{v,[u,w])} p 



{g{VitV - VyU, w) + g{v, V„w - Vt„u)} p 

where V„ is the Riemannian covariant derivative on the manifold {M,g). 
Thus 

G{uv, w) = y {V„ [g{v, w)]} p - J v^w^ [V,;Uj + VjU,;] p 

After an integration by parts the first term is zero, when u has zero diver- 
gence. It follows that 

[uv\j = — [ViMj + Vjlti] + Vicftiu, v) 

where (j){u, v)is to be chosen such that the Ihs has zero divergence: 

Thus, absorbing V(5(u, i;))into the gradient term, 

[uv + vu] = — [VtjTi + V„w] + \7ip{u, v) 



But 



so that 



[Duv] = ^{[u,v]+WyU + Wuv} + Wp(u, v) 



[u,v] = VuV - V„u 



[Duv] — \/yU + \7p[u, v), V^p(w, v) + divV^M = 0. 
In particular 

DyV = V„w + Vp 

which p chosen such that the divergence is zero. So the geodesic equation 
on the group of volume preserving diffeomorphisms is the Euler equation, as 
promised. 

The L^— inner product allows us to split the space of vector fields into the 
subspace of divergence-free vector fields (which is also a subalgebra) and another 
subspace of gradients: 
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u ^ 11^ + V(j){u), divu"^ = 0, V^(/)(u) = divu. 



This is an orthogonal decomposition. The Riemannian covariant derivative in 
the group of volume preserving diffeomorphisms is just the transverse projection 
of the covariant derivative on {M,g): 



Thus we can split the L^— inner product on vector fields into transverse and 
longitudinal pieces 



T{u, v) = G(u^, w^), S{u, v) = G(V0(w), V<^(w)). 
In other words, even if w is not of zero divergence 

G{DuV, w) {VuV, w) 

7.3 Curvature of the DifFeomorphism Group 

R{u, v) ^ G v) - G {DuDyU, v) + G {D.^DuU, v) 

= G v) +G {DyU, Duv) - G [DuU, Dyv) 

= G {y[u,v]U, v) + T (V-„w, V„u) - T (VuU, Vyv) 

where we use the fact that v has zero divergence. 
Now, 



R{u, v) ^ r{u, v)p + G (V„Vi,M - V^u V„u, v)+T (V^u, V„t;) - T (V„m, V„w) 



Set w = Vt,u, which may not have zero divergence, even though u and v are 
of zero divergence. Then, 



DuV 



G{u, v) = T{u, v) + S{u, v) 



where 



r(u, v)u + V„Vt,u — V„VijU 
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= -G{w, V„u) = -S{w, S/uv) - T{w, V„w) 
And similarly for G(Vt, V„u, w). Thus 

v) = r(u, v) + S'(V„u, Vuw) — 5'(V^u, V„w) 

where 




is the curvature form of {M,g) averaged by the density. It is not difficult 
now to work out explicit answers for a flat torus and recover Arnold's original 
results. We flnd the present form more useful as well as more general. 

8 Conclusion 

Arnold's geometric explanation of the instabilities of an ideal fluid raises an 
important physical question. Do the Euler equations describe the time evolution 
of real fluids? Any theoretical description of a physical system is an idealization 
in which small forces are ignored. This is usually justified because small forces 
have small effects. 

Near an unstable point this breaks down . Just think of a pendulum balanced 
precariously on its head, at the unstable equilibrium point. A small force in 
either direction can topple it in one direction or the other . The final outcome 
is random, since no one can predict with certainty such tiny perturbations. 

What if a system is dynamically unstable at every point in its phase space? 
This is the case with an ideal fiuid. In this case no where is it justified to ignore 
the unpredicatable small forces. Instead of predicting these forces, we must 
model them as a stochastic process; for example, as Gaussian white noise. The 
time evolution of the system in then a stochastic differential equation. 

What is the evolution of an unstable dynamical system under infinitesimally 
small white noise random forces? In a later paper [lO] we will show that this 
can be reduced to a deterministic dynamical system, but with twice as many 
degrees of freedom. The Euler equations will then get replaced with equations 
for geodesies on the tangent bundle of the diffeomorphism group. 
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